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Major Drivers of the Current Work 


• TOOLS FOR \FROSP \C\[ PhSIGN 

Decrease design cycle time => Rapid turn-around 

Increase design/process fidelity => High accuracy and low variation 

Increase discipline integration => Increased range of options via IT 


Turbo-pump component analysis ==> Entire turbo pump simulation 

• HPCC CAS Level I milestone : 

To demonstrate a 1000 times speed up in June 2001 over what was 
possible in FY92. 



Objectives 

• To enhance incompiv 

ssible How simulation capability h»r 

developing aerospace 

vehicle components, especially, unready 

flow phenomena associated . 

\ith high ^peed turbo pump 


199+. WOO 
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Shuttle Upgrade SSME-rig 1 
(code validation) 




Wet Guide Vane (IGV) 

•15 Blades 

•Pitch, p s 24 degrees 

•Blade Inlet Angle (mean), P^y, 3 90 degrees 
•Blade Exit Angie (mean), P, GVJ = -*5 degrees 

Clearance between KJV and Impeller, i = 0.12 inches 

InqjeOer 

•6+6+ 12 Unshrouded Design 
•Pitch, p a 60 degrees 

•Blade Inlet .Angle (mean), B^ r , » 23 degrees 
•Blade Exit Angle (mean), = 65 degrees 
•Clearance between blade LE and Shroud, r = 6.0056 Inches 
•Clearance between Made TE and Shroud, x = 0.0912 Inches 

Clearance between Impeller and Diffuser, r a 0.050 Inches 

DUDaer 
•23 Blades 

•PHch, p = 15.652 degrees 

•Blade Inlet Angle (mean), P^, = 12 degrees 

•Blade Exit Angle (mean), p^ = 43 degrees 


H 


' 4 / 




grid for tip 
clearance 


6 long blades / 6 medium blades /1 2 short blades hub 

60 Zones / 19 Million Grid Points grid 

Smallest zone : 74,025 /Largest zone : 899,248 

Overset connectivity is obtained by using DCF module OVERFLOW-P 


Page 5 







Page 6 




Time-Accurate Formulation 


• Time-integration scheme 

Artificial Compressibility Formulation 

• Introduce a pseudo-time level and artificial compressibility 

• Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 


Pressure Projection Method 

• Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 


Pressure Projection Method 


• Evaluate time-accurate features of two primitive variable methods 
designed for 3-D applications 

• Solve for the auxilary velocity field, using implicit predictor step 

— («; Vp n +h(uf) 

At 

• The velocity field at time level (n+1 ) is obtained by using a correction step, 

2 

— (u?"-u;)=-Vp’"'+h(ur')-Vp"+h(u") 

At 

• The incompressibility condition is enforced by using a Poisson equation for 
pressure ( P =p”*'-p n ) 

V 2 /?-— Vi i* 

At 
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Pressure Projection Method(INS3D-FS) 


• The discretization or the mass conservation equation in tinne volume 
form u la J ion > 

• New dependent variables, s »» 


u f = s* u 
u n = s' 7 " 

U ; = S u 

• Computing time : 

80 p-secs/grid point/iteration 

• Memory usage: 

70 times number of grid points 
in words 


r’/S.fc-O*. M/2 








1+1/2. K-V2 t*W 


Artificial Compressibility Method (INS3D-UP) 


Time-Accurate Formulation 

• Discretize the time term in momentum equations using second-order three- 
point backward-difference formula 

/ yi + l 

3U_ + dV_ + dW =Q 3q--4q- 

3£ 3?7 3£ 2A/ 

• Introduce a pseudo-time level and artificial compressibility, 

• Iterate the equations in pseudo-time for each time step until incompressibility 
condition is satisfied. 

^ ^pn + l.m+\ — pn+l.m^zz — qn + \,m+\ 

At 

L5 ( ..... i — 3 q—-4q-+q-‘ 

r 2A t 

• Computing time : 50-120 p-secs/grid point/iteration 

• Memory usage: Line-relaxation 45 times number of grid point in words 

GMRES-ILU(O) 220 times number of grid point in words 
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Impulsively Started Flat Plate at 90° 


• ft mo -r P'hih 

Artificial compressibility (INS3D-UP) : Line Relaxation 



Time (sec) 


Impulsively Started Flat Plate at 90° 


• Time of Stagnation Point 

Artificial compressibility (INS3D-UP) : GMRES 



1 1 j r 

0 2 4 6 

Time (sec) 
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Summary 


• An efficient solution procedure for time-accurate solutions of 
Incompressible Navier-Stokes equation is obtained. 

- Artificial compressibility method requires a fast convergence scheme.. 

- Pressure projection method is efficient when small time-step is required. 

- The number of sub-iteration is reduced significantly when Poisson solver is 
employed with the continuity equation. 

- Both computing time and memory usage are reduced (at least 3 times). 

- DCF module in OVERFLOW-D is incorporated with INS3D. 

- MP! /Open MP hybrid parallel code has been completed and benchmarked. 

• Work currently underway 

- Multi Level Parallelism (MLP) of INS3D. 

- Overset connectivity for the validation case (SSME-rigl) 

- Experimental measurements at NASA-MSFC. 

- Computational model for boost pump 
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